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In a topological insulator (TI) the character of electron 
transport varies from insulating in the interior of the ma- 
terial to metallic near its surface. Unlike, however, ordi- 
nary metals, conducting surface states in TIs are topo- 
logically protected and characterized by spin helicity 
whereby the direction of the electron spin is locked to 
the momentum direction. In this paper we review se- 
lected topics regarding recent theoretical and experimen- 
tal work on electron transport and related phenomena in 
two-dimensional (2D) and three-dimensional (3D) TIs. 



The review provides a focused introductory discussion 
of the quantum spin Hall effect in HgTe quantum wells 
as well as transport properties of 3DTIs such as sur- 
face weak antilocalization, the half-integer quantum Hall 
effect, s + p-wave induced superconductivity, super- 
conducting Klein tunneling, topological Andreev bound 
states and related Majorana midgap states. These prop- 
erties of TIs are of practical interest, guiding the search 
for the routes towards topological spin electronics. 
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1 Introduction. The situation when a material be- 
haves as a metal in terms of its electric conductivity and, 
at the same time, as an insulator in terms of its band struc- 
ture is extremely unconventional from the viewpoint of 
the standard classification of solids. That is why the re- 
cent discovery of a class of such materials - topological 
insulators (TIs) - has generated much interest (see e.g. re- 
views 0|2)). The dual properties of the TIs are especially 
well pronounced in two-dimensional (2D) systems which 
are also known as quantum spin-Hall insulators (QSHIs) 
|[3]|4][5]|6]|7]. In QSHIs, the metallic electric conduction is 
associated with propagating states that occur only near 
the sample edges, while the conduction in the interior is 
suppressed by a band gap like in ordinary band insulators. 
These edge states originate from intrinsic spin-orbit (SO) 
coupling and are profoundly different from those appear- 
ing in quantum Hall systems in a strong perpendicular 
magnetic field [8 9|. The key distinction lies in the role of 
the time reversal symmetry. In the QSHIs the SO coupling 
preserves the time-reversal symmetry, resulting in a pair of 
counter-propagating channels on the same edge as opposed 
to one-way directed (chiral) edge states in quantum Hall 
systems. Remarkably, the spin and momentum directions 
of the two QSHI edge channels are locked in the oppo- 
site ways so that these states are characterized by opposite 
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spin helicities and orthogonal to each other. As a result, 
such helical edge states have a nodal band dispersion (see 
also Fig. Q3 which is topologically protected against any 
structural or sample imperfections that do not cause spin 
scattering. 

The QSHIs have been experimentally realized in 
HgTe/CdTe quantum wells with inverted band structure 
and strong intrinsic SO splitting coming from the atomi- 
cally heavy mercury ||5]|6]|7). In these structures the helical 
edge states generate nonlocal transport effects in zero mag- 
netic field that have no analogues in the conventional 2D 
semiconductors, allowing one to detect the helical edge 
states in appropriately designed Hall bar devices |7). A 
variety of other properties of the QSHIs and, specifically, 
HgTe quantum wells have also got in the focus of re- 
cent theoretical [[T0][IT][n]^ 
|22l|23|24l|25]|26l|27l|28l|23 

|39]|40l|4D 42 43 44 45 46 47]|48l|49][50][5T]l52j[53l[54ll53 
I561I57 ,58 ,59 .60 ,6 Ill62l63l64l and experimental ||65ll66l 
I671l6^l6^l70ll7l1l72ll73ll74l research. 

The higher dimensional analogues of the QSHIs are 
the three-dimensional (3D) TIs ll731l76ll77ll7^l79ll80l[8T1 . 
In 3DTIs the topologically protected electronic states ap- 
pear on the surface of a bulk material. These surface states 
have a nodal band dispersion in the form of a Dirac-like 
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Figure 1 Schematic of band dispersion in 2D HgTe quan- 
tum wells: (a) ordinary band insulator with normal gap 
Ai > versus (b) topological quantum spin-Hall insu- 
lator (QSHI) with inverted gap Ai < 0. In the latter case, 
two crossing spectral branches correspond to a pair of edge 
states with opposite spin helicities r = ±1. We use Hamil- 
tonian © with A = 380 meV-nm, \M\ = 10 meV, 
B = T> = and boundary condition ©. 



Figure 2 Energy bands (in meV) of a surface state in a 
3DTI versus in-plane wave-numbers k x and k y (in nm -1 ) 
from effective Hamiltonian (f34t . The Fermi level lies in the 
conduction band at E = 0. We chose A = 300 meV-nm, 
W = 50 meV-nm 3 and T> = meV-nm 2 (adapted from 

ESI). 



cone reflecting a continuum of momentum directions on 
the surface (see also Fig. |2}. The family of materials and 
heterostructures which can host surface states is pretty 
large. They were first predicted in inverted semiconduc- 
tor contacts [82,83 1. The coexistence of the metallic sur- 
face states with the bulk gapped band structure, i.e. the 
3DTI phase, has been established theoretically for the 
semiconducting alloy Bii^Sbz ll76l . strained 3D layers 
of a-Sn and HgTe IT761 . the tetradymite semiconductors 
Bi2Se3, Bi 2 Te3, and Sb2Te3 1811 . thallium-based ternary 
chalcogenides TlBiTe 2 and TlBiSe 2 I84ll85ll86l as well as 
Pb-based layered chalcogenides 1 87 , 88 1 . Experimentally, 
topological surface states have been observed by means 
of angle-resolved photo-emission spectroscopy (ARPES) 
in Bh- x Sb x M, Bi 2 Se 3 (83, Bi 2 Te 3 flU, TlBiSe 2 
EH52I531, TlBiTe 2 J93), PMBi^Sb,)^ HI and 
PbBi 2 Te 4 E3. 

Beside the ARPES and band structure calculations, 
there has been a growing number of experiments (96 , 97 , 
l9lll9 W00l[l^[l^[T0^T04ll 1 0511 1 06lfT07ll 1 0811 1 0911 1 1 01 
fTmil 12111 13111 14U115U116l[TT7lll 18111 19U120lll21lfl22l 
and theoretical studies 1 1 2311 1 2411 1 2511 1 2611 1 27lfT28lfl29l 
1 1 3011 131 irT32ll 1 3311 1 34U 1 3511 1 3611 1 371fT38lfT39ll 14011 14T1 
1 1 4211 1 43lfT44ll 1 45 II 1 4611 1 47ll 1 4811 1 49lfl30lfT5T1[T52l[T53l 
ll541ll55irT56Tll57U158U159lll60lll61irT62lfT63lll641ll65l 
166 167 168, T69ll7O0171lll72M173lfr74lfT75lll76M177l 
[178111 7911 180lfT8TH devoted to helical transport in 3DTI 
materials. Like in the QSHIs, the surface charge carriers 
in 3DTIs are characterized by a well-defined spin helic- 
ity, i.e. the locking of the spin and momentum directions. 
This has several implications for electron transport which 
will be the subject of the present review. First, the helicity 



conservation along closed electron trajectories prevents 
ordinary electron localization on the surface. Instead, ex- 
periments 19611 1 0311 1 05lfT04ll 1 1 1 II 1 06lfT07TfT08l[TT9l[T20l 
and theoretical calculations 0155lll56lfT57lll58lll78IIT79l 
have indicated that the surface quantum transport in the 
3DTIs should exhibit weak antilocalization with a positive 
magnetoresistance. In thin TI films with surface and bulk 
states, both weak antilocaliation and localization regimes 
have been predicted 1158111791 . Second, in a strong per- 
pendicular magnetic field the unconventional half-integer 
quantum Hall plateaus are expected for a single spin- 
helical surface ll7^[n4l[T34l[T37l[T62l . as has been found 
in graphene 01821 and in zero-gap HgTe/CdTe quantum 
wells 0681 . Experimentally, unusual odd and even quantum 
Hall plateaus have been observed in epitaxially strained 
3DTI HgTe 11 101 . which may indicate the contribution of 
two spin-helical surfaces, one at the top and one at bottom 
of the sample. Beside the magnetotransport, the quantum 
Hall dynamics in 3DTIs manifests itself in the Faraday 
effect at THz frequencies 11221 . Another interesting impli- 
cation of the spin helicity is the possibility of unconven- 
tional surface superconductivity. It is expected to occur in 
proximity of a singlet s-wave superconductor (e.g. Nb or 
Al) which induces both singlet s-wave and triplet p-wave 
correlations, as a result of the broken spin-rotation sym- 
metry 1123lll30lfT32lll49lll50l[T70l . As a signature of the 
p-wave correlations, topological Majorana midgap states 
have been predicted in superconductor (S)/3DTI junctions 
(see e.g. [123,2, 169, 170|). The theoretical search for un- 
conventional superconductivity in 3DTIs goes in parallel 
with experimental progress in fabrication and characteri- 
zation of such hybrid superconducting systems 11 13111 14l 
II 15111 16l fTT7lfT2Tl . 
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Figure 3 Band structure of a HgTeZHgo.3Cdo.7Te quan- 
tum well from k • p model: (a) electron El, E2, ... and 
heavy-hole HI, H2, ... subband energies versus well thick- 
ness d, (b) in plane dispersion at the critical thickness 
d c w 6.3nm, and (c) a 3D plot of the Dirac-like low-energy 
spectrum for d — d c near the r point of the Brillouin zone 
(adapted from l68l ). 



As transport in TIs is a multifaceted problem, this re- 
view aims to discuss it from different angles, starting from 
the QSH effect and related phenomena in 2D HgTe quan- 
tum wells and continuing to the surface weak antilocal- 
ization, half-integer quantum Hall effect and s + p-wave 
induced superconductivity in 3DTIs. Regarding the super- 
conducting properties, we introduce a generalized Fu-Kane 
model correctly accounting for the energy-dependence of 
the proximity effect in 3DTIs. We then analyze the su- 
perconducting Klein tunneling, topological Andreev bound 
states and related Majorana midgap states as well as the 
fractional AC Josephson effect 111 8311 184fl in such proxim- 
ity S/TI/S junctions. 

2 Topological insulators in two dimensions: 
Quantum spin Hall insulators. 

2.1 Effective Hamiltonian for HgTe QWs. We will 
discuss the QSHIs in the context of their experimental 
realization in HgTe/CdTe quantum wells (QWs) (SHU- 
HgTe is a zinc -blende-type semiconductor whose band 
structure can be described by the eight-band Kane model. 
Unlike conventional zinc-blende semiconductors, due to 
large relativistic corrections HgTe has an inverted band 
structure where the Jg band (originating from metallic s- 
orbitals and usually acting as the conduction band) has a 
lower energy than the Jg band which derives from chalco- 
genide p-orbitals 01851 . Consequently, in HgTe QWs with 
large enough thickness d the energy subbands are also in- 
verted (see Fig.[3k). In Fig. [3] the Ig-derived subbands are 
denoted by HI, H2, ... because they correspond to heavy- 
hole subbands in non-inverted semiconductors, whereas 
the electron-like Jg -derived subbands are denoted by 
E1,E2, .... With the decreasing QW thickness d, the en- 
ergies of the E subbands increase as a result of quantum 
confinement, whereas those of the H subbands decrease, 



as shown in Fig. [3^. This results in the normal band se- 
quence in narrow QWs. The different d-dependences of 
El and HI subbands imply a critical thickness, d c , at 
which the band gap is closed. For d « d c and near the 
r point, an effective four-band model involving double 
(Kramers) degenerate El and HI subbands can be de- 
rived from the eight-band Kane model l4 l26l . Introducing 
basis states \E1+), \H1+), \E1-) and \H1-) (where 
± denotes Kramers partners), one can write the effective 
four-band Hamiltonian for the system as follows I4l l26l : 



H 



HgTe 



h(k) 
h*(-k) 



h(k) = A{a x k x - (Jyky) + MkCTz 
Mk = M+Bk 2 . 



(1) 



Vk 2 a , (2) 
(3) 



The two diagonal blocks of Hn g Te O describe pairs of 
states related to each other by time reversal symmetry 
(Kramers partners). Each of the blocks has a matrix 2x2 
structure with Pauli matrices u x . v . z and unit matrix do rep- 
resenting the two lowest-energy subbands El and HI. The 
linear terms in Eq. (fJI (proportional to constant A and in- 
plane wave-vectors k Xt y) describe the El-Hl hybridiza- 
tion, while A^k yields the band gap M at the r (k = 0) 
point of the Brillouin zone. The positive quadratic terms 
Bk 2 and T>k 2 take into account the details of the band 
curvature in HgTe QWs |4i]. The Hamiltonian (Q~|i can be 
extended to include the spin-orbit coupling between the 
Kramers partners Il26l l6ll . 

Using unitary transformation H — > UHW with U = 

( -ia 'tf ) > we can cast tne Hamiltonian ([TJ into a Dirac- 
like form 



H 



[Ak + M k z) + £>k 2 T cr , 



(4) 



where Pauli matrix r 2 and unit matrix To act on the 
Kramers partners. We note that despite the effective mass 
term A4ut z ct z the Hamiltonian (0]i is invariant under time 
reversal, i.e. HT = H, where T — ir y a x C is the time- 
reversal operator, with C denoting complex conjugation. 

2.2 Simple analytic model of the QSH insulator. 
The QSHI state is realized when a QW with an inverted 
gap Ai < is sided by ordinary band insulators with 
Ai > (see also Fig. QJ)). The system boundaries play 
a role of topological defects - "mass" domain walls - that 
bind electronic states near the edge so that they decay on 
both sides of the boundary. In this subsection we discuss 
specific properties of such edge states: 

-the edge-state spectrum is gapless and merges into the 
bulk spectrum above the band gap; 

-the edge-states are the orthogonal eigenstates of the he- 
licity operator E = t z ct ■ k, where k is the unit vector 
in the direction of the edge-state momentum. For this 
reason, the QSH edge states are called helical; 
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-local static perturbation V preserving time-reversal 
symmetry does not couple the QSH edge states. 

In order to illustrate these properties we will make two 
simplifications. First, we will omit all the terms cx k 2 in 
Hamiltonian ©. This is justified since the edge states oc- 
cur in the vicinity of the _T(k = 0) point. Hamiltonian (|4]i 
takes, then, the form H = t z <j ■ (.Ak + Adz), which in po- 
sition representation corresponds to the following equation 
for the four-component wave function ^(r): 



[ea - T z a ■ {-iAV + Mz)]&(r) = 0. 



(5) 



Second, we will assume that our system is confined by a 
normal band insulator with the infinite mass, M — s- +oo. 
It is known [1861 that such infinite mass confinement can 
be modeled by an effective local boundary condition which 
in our geometry (e.g. at y = 0) reads 



${x, y = 0) = T Q a x $(x, y = 0). 



(6) 



This boundary condition is specific to Dirac fermions with 
linear spectrum. It ensures vanishing of the normal compo- 
nent of the particle current without putting ^(x, y) to zero 
at the boundary. The use of the infinite mass confinement 
(O is complementary to the tight-binding calculations (see 
e.g. Refs.|6j .18]) and other continuum models of the heli- 
cal edge states, which include the quadratic (k 2 ) terms and 
employ the hard wall boundary conditions (see e.g. Refs. 

CUE]). 

We seek solutions to Eq. (0 in form of the two eigen- 
states, <Pk ,± ( r )> °f diagonal matrix t z uq propagating along 
the edge (in the x-direction) and decaying exponentially 
away from it (in the y-direction): 



(7) 
(8) 



with a real positive decay length A > 0. The symbol ® de- 
notes a tensor product of an eigenstate of t z (first column) 
and the wave function in a space (second column). The 
conditions for the nontrivial solutions for the coefficients 
x b\k± and x S>2k± follow from Eqs. (0 and ©, yielding two 
equations for A and e: 



A- 2 -A/ 2 A4 2 = fc 2 -e 2 A4 2 , 
X- 1 + AI/A = k - e/Ar, t ■■ 



±1. 



(9) 
(10) 



We notice that the left-hand-side of Eq. (TTOb does not con- 
tain index r, whereas the right-hand-side does. This can 
only be true if both sides of Eq. ( TTOb [and those Eq. (0] 
vanish independently, which yields a solution with a gap- 
less linear dispersion and a real decay length: 



e kT = Akr, A = -A/M, 



M < 0. 



(11) 



Since A must be positive, the edge states exists only in a 
system with the inverted negative gap, disappearing when 
A4 turns positive. Their propagation velocity v = A/K co- 
incides with that of the bulk states above the gap (see also 
Fig. [TJ. In a narrow QSHI the overlap of the edge states 
from the opposite sides results in a gapped edge-state dis- 
persion IfTTI . 

The edge-state wave functions normalized to half- 
space < y < oo are given by 



*Jb,+ (r) = 



#*,-(') = 




ikx— \M\y/A 



ikx—\M\y/A 



(12) 



(13) 



The key feature of the edge states (Tl2l) and ( Tl3b is that 
they are orthogonal eigenstates of the helicity operator 

E = t z u x : 



E& kiT (r) =T!P fc , T (r), 



r = ±1. 



(14) 



The helicity S is defined as the projection of vector £ = 
T z cr on the direction of the edge-state momentum k\\x. 
Since the matrix structure of E derives from the SO-split 
energy bands, it is also called the spin helicity. Equation 
(fT4l is a manifestation of the time-reversal symmetry and 
the fact that the QSH state is generally characterized by a 
Z2 topological invariant (3J, It is easy to see that one helical 
channel can be obtained from the other by simply applying 
the time-reversal operator T = ir y a x C: 



TflU 



k.+ - 



(15) 



The helicity [see Eq. (TPfl il protects the edge states from 
local perturbations that does not break the time-reversal 
symmetry. Concretely, let us consider a perturbation V 
which is (up to a phase) invariant under time reversal: 



T t VT = V*. 



(16) 



Using Eqs. (TOT l and ( TToT l we can transform the matrix ele- 
ment (—k, —\V\k, +) as follows 

(-k,-\V\k,+) = -(k,+\T f VT\-k,-) = (17) 

= -<fc,+|y|-fc,-) = -<-fc,-|yt|fc,+>, (18) 

i.e. for a hermitian V = its matrix element is zero: 



(-k,-\V\k,+)=0. 



(19) 



In particular, spin-independent potential disorder cannot 
cause scattering between the helical edge states and their 
localization. The physics of localization in the QSHIs has 
been studied by using both field-theoretical [fJO 24| and 
numerical (see e.g. lfT8l l35,37,52,5H) techniques. In par- 
ticular, Ref. |[l"8"l has found that the edge backscattering and 
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a) QH edge transport b) QSH edge transport 





Figure 4 Schematic of edge transport in four-terminal de- 
vice: (a) chiral edge states in a quantum Hall (QH) sys- 
tem versus (b) helical edge states in the quantum spin Hall 
(QSH) insulator. 



where Tji is the transmission probability from contact i 
to contact j. For a chiral QH edge channel, Tji connects 
the neighboring contacts only in one propagation direction 
(see also Fig.|4^), such that 



T(QH) 



i+l,i 



1. 



i = l,...,N, 



(21) 



where N is the number of the terminals (e.g. TV = 4 in 
Fig.|4]i, with the convention that Ijv+i,jv = Ti.iv describes 
the transmission from terminal N to terminal 1 . Assuming, 
for concreteness, that the current flows from terminal 1 to 
terminal 4, while leads 2 and 3 are used as voltage probes, 
we have 



magnetoresistance can occur as a combined effect of bulk- 
inversion asymmetry, sufficiently strong potential disorder 
and an external magnetic field. 

Another question regarding the role of the time- 
reversal symmetry is what happens to the QSHIs in a 
strong quantizing magnetic field? In this case the helical 
edge states appear within the gap between the highest hole 
(— |M|,r = +1) and lowest particle (| M |, r = -1) Lan- 
dau levels JT9|. The edge-state dispersion is still nodal, but 
nonlinear, i.e. the two crossing spectral branches have now 
different group velocities. The latter circumstance make 
the helical edge states prone to get coupled by spin-flip 
scattering, which also generates the edge magnetoresis- 
tance Ifl9ll50l . 

We should emphasize that Eqs. ( fT7T > - $1~9[ do not gen- 
erally hold for interacting helical channels where the elec- 
tron backscattering may occur provided that the axial spin 
rotation symmetry of Hamiltonian (Q~|) is broken, e.g. by 
the Bychkov-Rashba spin-orbit coupling 1 20, 47, 60). The 
spatially random Bychkov-Rashba coupling has been pre- 
dicted to cause localization of the edge states in the pres- 
ence of weakly screened electron-electron interaction l20l . 
Ref. l)60l has shown that weak electron-electron interac- 
tion in the absence of the axial spin rotation symmetry al- 
lows for inelastic backscattering of a single electron, ac- 
companied by forward scattering of another. The inelas- 
tic phonon-induced backscattering in the helical liquid has 
been considered in Ref. l47l . 



2.3 Helical versus chiral edge transport. Nonlo- 
cal detection of the QSHI state. The helical edge trans- 
port in a QSHI differs markedly from the chiral edge trans- 
port in QH systems. In order to illustrate this we follow 
Ref. [7i| and consider a four-terminal device shown in Fig. 
|U Using the Landauer-Biitikker approach, we express the 
current, injected through contact i in terms of voltages 
Vj induced on all contacts as follows: 



TV 



.^(TjiVi-TijVj), 



(20) 



h = hi = -r(Vi - V 4 ), V 2 = V u V 3 = V 2 , h = -h- 
h 

(22) 

This yields a finite two-terminal resistance i?i4,i4 = 
Vl Y^ = h/e 2 and zero four-terminal (nonlocal) resis- 
tances R u i 2 = JSpVa = o, R u 23 = = and 

-*14 J 14 

i?14,13 = ^ = 0. 

In contrast, in a QSHI the helical edge channels con- 
nect the neighboring contacts in both propagation direc- 
tions (see also Fig.|4j)), such that 



T(QSH) i+hi = T(QSH)i, i+1 = 1, i = 1, N, 

(23) 

with the conventions Tjy+i jv = 7i,iv and TV.jv+i = 
Tjy l- Consequently, for a current flowing from 1 to 4, we 
find 

e 2 

h =h4 = j (2Vi -Vi- V 2 ), h = -h , (24) 
2V 2 - V x - V 3 = 0, 2V 3 -V 2 -V 4 = 7 (25) 

which yields the two-terminal resistance 

Vi-Va 3 h 



R 



14,14 



(26) 



hi 4 e 2 
and the four-terminal resistances Q 



Vi - V 2 _ V2-V3 1 h 

-^14,12 = Z = -Rl4,23 = j — (27) 

hi hi 4e^ 



R 



14,13 



V1-V3 
hi 



'14 

1 h 

2^' 



(28) 



The nonzero non-local resistances d27b and d28l ) are unique 
to the QSHI state, allowing its unambiguous experimental 
detection (7). It should be emphasized that the universality 
of the non-local resistances d27l i and d28l ) is just the con- 
sequence of the time-reversal symmetry and, therefore, is 
also expected for other proposed realizations of the QSHIs, 
e.g. in inverted InAs/GaSb quantum wells ll58ll . 
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Figure 5 Backscattering rate rT^ [in units of elastic sca- 
tering rate r _1 , see Eq. d30l ll versus carrier density n; 
A = 380 meV-nm and B = 850 meV-nm 2 (adapted from 



155|). 



2.4 Helical carriers and weak antilocalization in 
n-type HgTe quantum wells. The absence of the edge 
localization in the QSHI regime (i.e. when the Fermi level 
lies in the band gap) is closely related to the weak antilo- 
calization (WAL) effect observed when the Fermi level is 
pushed above the gap into the bulk conduction band. The 
latter case corresponds to n-type HgTe quantum wells in 
which charge carriers behave as a 2D helical metal. Assum- 
ing that it is described by the same Hamiltonian (@}, we can 
calculate the disorder-induced quantum-interference cor- 
rection 5a to the classical Drude conductivity (for more 
details see Sec. [3]). In the leading logarithmic order Sa is 
given by 115511 



2e 2 t~ 
8a xx (n,M) f« — In—, — r 

T kA +Tu; 



! 2 / 7s - 7T 



' .VI 



M ~ '<P 

2 



JB = 7T 1 



M 



Bkp 



\jA 2 k 2 F + (M + Bkp) 2 



(29) 

(30) 

sjlixn. 
(31) 

where r is the elastic scattering time and t v is the de- 
phasing time. The third new time-scale t_m appears due 
to the lack of topological protection of the gapped spec- 
trum: its geometrical Berry phase 73 OTT ) no longer coin- 
cides with the universal value ir characteristic of a massless 
Dirac cone 111871 . The presence of the gap Ai + Bkp en- 
ables scattering between opposite-momenta states k^ and 
— kp (backscattering) on the 2D Fermi surface (see also 
Ref. [30]). Thus, r^ 1 is the rate of such backscattering. 
It is an interesting band-structure effect (see also Fig. [5]): 
for the normal band structure with M > the rate tJ^ 
is finite at any carrier density n, whereas for the inverted 
band structure with M. < the backscattering rate t^ 1 
is strongly nonmonotonic with a zero at a specific carrier 
density n = —M/2nB. 

The role of the geometrical phases in transport 
quantum-interference effects in HgTe quantum wells has 



recently been studied numerically in Ref. Il62ll . In Ref. 
l63l analytical results for the WAL in HgTe quantum wells 
have been obtained with the account of both bulk-inversion 
asymmetry and Rashba SO coupling. Experimentally, the 
WAL in HgTe has been observed in 2D quantum wells 
[66 74] and in strained 3D layers 111 1 Q . Also, experiments 
have revealed cyclotron resonance phenomena for the he- 
lical carriers in HgTe quantum wells 169117 II . 

3 Spin-helical transport in 3D topological insu- 
lators. The ID edge states of HgTe/CdTe quantum wells 
discussed above have higher-dimensional analogues, such 
as gapless states on a 2D interface between two bulk semi- 
conductors with normal and inverted gaps 18211831 . Gen- 
erally, a 3D system with gapless surface states occuring 
inside the bulk band gap realizes a topogically nontrivial 
insulating state, the 3DTI l75ll76ll77ll78l[8l1[Tir2l . Since the 
gapless states in 3DTIs appear only on the surface, they 
are exempt from the fermion doubling theorem [188]. For 
this reason, the topological surface state consists of an odd 
number of Weyl-like fermions, each described by a two- 
component spinor wave function. In Bi2Se3 and Bi2Te3 
there is a single Weyl-like fermion species I8l l89ll90l ) 
which can be described by an effective 2D Hamiltonian 
I189019OI : 



H = A(a x k y 



CTyk x ) 



it- 



(32) 



(33) 



where k = (k x , k y ) is the wave vector on the surface and 
A, T> and W are band structure parameters. Upon unitary 
transformation H -> UH W with U = (01), Hamilto- 
nian (l32l takes the form: 



H = a ■ (Ak + W k z) + Vk 2 (T , 



(34) 



which looks similar to the HgTe quantum well Hamilto- 
nian There are, however, two important distinctions 
between these Hamiltonians. First, here the basis functions 
correspond to | and — ^ electron spin projections, i.e. Pauli 
(tj x ,y,z) an d unit (do) matrices act on real spin indices. Sec- 
ond, the yVg.-term in Eq. (l34l l is cubic (odd) in momen- 
tum k, causing no gap at k = 0. This term does not break 
k, <r — >• k, —a invariance, which is the real time-reversal 
symmetry in this case. Instead, it causes hexagonal warp- 
ing 1191111921 of the surface spectrum (see, also Figs. [2]): 



E(k,<p a ) = Vk 2 ±\ A 2 k 2 + 



W 2 k 6 



(l + cos60 n ), (35) 



where the angle <p n indicates the momentum direction. 

3.1 Persistence of surface spin-momentum lock- 
ing in the presence of disorder. In Eqs. d32j and (l34t 
the first linear term is analogous to the Bychkov-Rashba 
spin-orbit interaction that occurs in a conventional 2D elec- 
tronic system confined in an asymmetric quantum well 
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111931 . The key difference here is the strength of the cou- 
pling constant A. For the topological surface states, the ex- 
perimentally determined coupling constants typically are 
A = 355 meVnm forBi 2 Se 3 |fT92l and A = 255 meVnm 
for Bi2Te3 f90l l9U. These values are really giant com- 
pared to the Bychkov-Rashba coupling o.br ~ 1 -j- 10 
meV nm in conventional GaAs-based 2D electronic sys- 



tems (see e.g. Ref. 119411 ). This large quantitative differ- 
ence has two important physical consequences. The first, 
as mentioned above, is the formation of the conical band 
dispersion [see also Fig. [2). The second is the robustness 
of the surface spin-momentum locking against potential 
impurity scattering (and other spin-independent scattering) 
as long as the energy separation between the valence and 
conduction bands [see Eq. < f33T >1 is larger than the spectrum 
broadening Ti/tq: 



2jA 2 k 2 F 



= 2(E F -Vk 2 F )^>h/T , (36) 



where tq is the elastic life-time, kp and Ep is the Fermi 
momentum and energy (the Fermi level is, for concrete- 
ness, located in the conduction band). Due to the large con- 
stant A the condition (f36b can be met already for the Fermi 
energy values, Ep, of several tens of meV. Moreover, the 
condition ( T36l l simultaneously implies a diffusive metallic 
transport regime in which 



2Ep To/h w kpvpT = kpl 3> 1, 



(37) 



because for a conical dispersion the transport momentum 
relaxation time r w 2to 111871 (v F is the Fermi velocity 
and i = vpT is the transport mean-free path). Thus, the 
large values of the coupling constant A [see Eq. ( T36l ll en- 
able the diffusive metallic regime d37l i without significant 
spin relaxation, i.e. the spin-helicity of the carriers is pre- 
served during their diffusion on a disordered TI surface. 
For Hamiltonian d34l i the conservation of the spin-helicity 
can be expressed as 



,4k + >V k z 
|.4k + W k z| 



const, 



(38) 



where the constant h takes the values ±1 for the con- 
duction and valence bands, respectively. The absence of 
the spin relaxation under condition d3~8l defines a distinct 
transport regime - the 2D helical metal -, opening hitherto 
unexplored routes for theoretical investigations of disor- 
dered systems (e.g. Refs. E41I140U152II163II172II ). 

3.2 Surface weak antilocalization in 3D TIs. One 
important example of observable transport phenomena in 
disordered systems is the weak localization (WL) effect 
II 1 9511 196ll 1971 . It is associated with closed classical elec- 
tron trajectories that self intersect with opposite electron 
momenta, k and — k (see Fig.|6h). On such trajectories, two 
states which are related to each other by time reversal in- 
terfere constructively, yielding, in the absence of the spin- 
orbit coupling, a negative WL correction Sa to the classical 




enclosed area 



Figure 6 (a) Schematic of closed electron trajectories giv- 
ing rise to weak localization [see also Eq. ( f39t and text]. 
Shaded areas indicate a space element of the trajectory, 
\FVF<lt, and the area enclosed by the trajectory, Dt, which 
we use to estimate the return probabilty dP re t{t) in Eq. 
(l39l . (b) Same trajectories on the surface of a 3D TI in- 
volve opposite momenta, k and — k, and opposite spins, er 
and — <t, as a result of the conservation of the spin-helicity 
[see also Eq. (HTI) and text]. This leads to the weak antilo- 
calization conductivity d42b . 



Drude conductivity ap>. The relative change 8a jap, is the 
measure of the classical return probability integrated from 
the shortest diffusion time-scale r to the largest one, given 
by the dephasing time r v : 

— ~ ~ r dP r^) ~ (39) 



WDty 

The negative sign here is the consequence of the fact that 
the particle returns to its original location with the opposite 
momentum k [i.e. the Kubo formula for the conductivity 
correction contains the product of the two opposite-sign 
velocities]. The return probability dP re t(t) is proportional 
to a space element of the trajectory, Xpvpdt, divided by 
the typical area {^/Dt) 2 enclosed by the particle trajectory 
at time t (D is the diffusion constant). The conductivity 
correction d39i > is logarithmically divergent and has a uni- 
versal prefactor 111 9511 196|[T97l : 



5a 



2e 2 , t v 

In — . 

2nh t 



(40) 



where the factor of 2 in the numerator is due to the spin 
degeneracy. 

In contrast to Eq. (|4D1 . according to experiments 1%1 
ll03lll05lfT^[TTT1[l^[l^[T08l[TT9in^ and theoretical 
calculations 0155lll56|[T57lll58lll78|[T79ll , the topological 
surface states in 3D TIs exhibit the weak antilocalization 
(WAL) effect characterized by a positive conductivity cor- 
rection Sa > 0. The change of the sign can be easily 
explained by the fact that the classical motion along the 
loop is now subject to the conservation of the spin-helicity 
[see Eq. (l38ll. Indeed, upon returning to its original loca- 
tion the particle with momentum — k must have the same 
helicity as it had initially with momentum k. This requires 
the change of the direction of the electron spin from a to 
— er (see also Fig. [6J?) so that we have the identity 



,4k + W k z 
|.4k + W k z| 



• ;^ + . (4i) 



\Ak + W k z| 
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We keep only the spin-singlet Cooperon, C(q) = l/(Dq 2 + 
t~ 1 ) , because this is the only gapless cooperon mode, with 
D = v f t/2 being the diffusion constant [we also intro- 
duce an auxiliary quantity Dq = {tq/t)D\. The factor 
(t/to) 2 originates from the independent renormalization 
of both current vertices (see Fig. [7}:) due to the difference 
between the elastic life-time To and transport relaxation 
time t for helical carriers 11 1551 : 



Figure 7 Diagrammatic representations for (a) bare and 
dressed Hikami boxes for the correction to Drude conduc- 
tivity, (b) Bethe-Salpeter equation for the Cooperon, and 
(c) equation for the renormalized current vertex in the lad- 
der approximation. Thick lines denote disorder-averaged 
Green's functions in self-consistent Born approximation, 
dashed lines - the potential (spin-independent) disorder 
correlation functions. 



The spin rotation cr 
Eq. <E3: 



-cr yields another "minus" sign in 



5a 



~~r In — . 

2-Kh T 



(42) 



This can also be viewed as the result of the it Berry phase 
accumulated along the trajectory loop [ 1 87 , 198|. 

The explicit calculation of 5a with the account of the 
warping Wt under condition d36b was carried out in Refs. 
II155II178I . We note that the conductivity correction 5a (l42l 
has the same form as for a conventional 2D electron system 



with spin-orbit impurity scattering 11991 or with Bychkov- 
Rashba and Dresselhaus spin-orbit interactions (see e.g. 
Refs. 1200112011 ). The reason is that the surface states of 
3DTIs with hexagonal warping and conventional 2D elec- 
tron systems with spin-orbit impurity scattering, Bychkov- 
Rashba or Dresselhaus spin-orbit interactions belong to the 
same - symplectic - universality class of disordered sys- 
tems. 

In contradiction with the aforementioned symmetry ar- 
gument, Ref. 1 160] has found a nonuniversal prefactor in 
Eq. ( 1421 which depends on the hexagonal warping strength. 
It is therefore instructive to demonstrate here that there is a 
cancellation of different warping terms 11551 . which yields 
a universal prefactor e 2 /2nh in Eq. (l42l . We will eval- 
uate the standard formulas for the quantum-interference 
conductivity correction, which are expressed diagrammat- 
ically in Fig. |7] (see e.g. Refs. |fl"97l , 112021 ). treating the 
warping < f33T > as weak perturbation onto the isotropic spec- 
trum under condition 



W 2 k%/2A 2 < 1. 



(43) 



Integrating out the single-particle momenta, we can ex- 
press the first diagram in Fig.|7^ - the bare Hikami box - in 
terms of the spin-singlet Cooperon C(q) as follows: 



da 



bare 



e 2 D 



Vt / 



dq 



C(q). 



(44) 



To 



1 



W*JA*k* F 



1 + W 2 k%/2A 2 ' 



(45) 



Here the squared warping term is averaged over the direc- 
tions of the momentum, k^, on the Fermi surface [(...) de- 
notes the averaging]. For helical carriers, dressed Hikami 
boxes (second and third diagrams in Fig. |7^) need also to 
be taken into account. Evaluating with the same accuracy 
with respect to W 2 k F /2A 2 <C 1 the second diagram in 
Fig-EK we obtain 



5a 



e 2 D fr_\ 21 



W 2 /A 2 k 



dq 



C(q). 
(46) 

In the absence of the warping, 5a dres equals to (—1/4) x 
fi a bare m a g reemen t with the calculations of the WL in 
graphene 1 202 1 . The third diagram in Fig.|7^ gives the same 

fi a dres as m £q ( |4g] | so t }j at me net resu lt j s 

5a bare + 25a dres 



5a 



e 2 D rr^l + W 2 F /A 2 k 2 F 

ttTi \Tq 



dq 



■C(q) (47) 



e 2 Dp 



TO 



(2^ C(q) = 2^ ln V 



3.3 Detecting surface states by WAL magneto- 
transport. So far we have treated the surface states in 3D 
TIs as purely two-dimensional. There is, however, a finite 
length A of order of a few nm over which they penetrate 
into the bulk of the material. Since A is quite small, the 
surface WAL conductivity is sensitive to the orientation of 
an external magnetic field with respect to the surface of the 
material, which could be used in practice to detect the sur- 
face states (see e.g. recent experiment [ 104|). Defining the 
magnetoconductivity as Aa(B) = 5a(B) — 5a(0) , one 
can obtain the following B-field dependences for perpen- 
dicular (_L) and parallel (||) field orientations: 



Aa±(B) 



e 

2nh 



B± = 



4|eK 2 ' 
(48) 



Aau(B) = -— In 1 
11 v ' 2-Kh \ 




B„ = 



V2\e\X£ 2 v 

(49) 

These equations contain the same phase-coherence length 
£ v = \/Dt v . However, the magnetic -field scales, B± 
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and B,, , on which the magnetoconductivity decreases, 
are distinctly different. The field B± corresponds to the 
Aharonov-Bohm magnetic flux of order of the quantum 
h/e through a typical area ~ S? enclosed by the interfering 



In order to identify the surface state one should extract 
its penetration length A from Eqs. (l48b and (l49b . Excluding 
i v from Eqs. d48b and d49l we can express A in terms of 
two parameters, B± and B., , which can be obtained inde- 
pendently from fitting the corresponding experimental data 
US): 



trajectories 11951 . while By corresponds to the same flux 
h/e, but through a significantly smaller area ~ \l v which 
is proportional to the surface-state penetration length in 
the bulk 12031. 



A = 




(50) 



The knowledge of this penetration length also helps to es- 
timate the critical thickness of the sample at which the two 
surface states start to overlap and the TI state disappears. 

To conclude this section, we should mention that the 
WAL in 3DTIs is suppressed and may even turn into the 
WL if an energy gap opens in the surface spectrum, as a re- 
sult of time-reversal symmetry breaking 11571 . This should 
apply to 3DTIs with magnetically doped surfaces (see e.g. 
Ref. 1 104]). The WL behavior has also been found in ultra- 
thin TI films in which the lowest bulk-state subbands are 
described by a massive 2D Dirac model [158]. Besides, 
quantum transport in electrically gated TI films may de- 
pend sensitively on the coupling between surface and bulk 
states 11791 . Such a coupling introduces new time scales on 
which the crossover from WAL to WL may occur 11791 . 

4 Quantum Hall effect on topological surfaces. 

Another manifestation of the surface spin-helicity in 3DTIs 
is the unconventional half-integer quantum Hall effect 
(QHE) (see e.g. Refs. 17^[n4l[ll4l[LT7lfT6a ). Such un- 
conventional QHE was first observed in single atomic lay- 
ers of carbon - graphene - 11821 . and has been regarded as 
a purely 2D phenomenon. The search for the half-integer 
Hall quantization in 3DTI systems is therefore a new chal- 
lenging task (see e.g. Refs. im fL34lll37U162|[TT0l ). Our 
discussion of the QHE in 3DTIs follows closely Refs. 
H137J162J ■ 

4.1 Landau quantization of surface states and 
conductivity tensor. Using a linear Hamiltonian H = 
vcr ■ (—ihV — eA) for a single topological surface state in 
a perpendicular magnetic field (A is the vector potential), 
we arrive at the following eigenvalue problem: 



H\n) = e\n), H 



~iV2hn B a 
iy/2hf2 B a) 



(51) 



where |.„) denotes an eigenspinor, a) and a are the rais- 
ing and lowering operators of a harmonic oscillator, re- 
spectively, and TiQb = %v/Ib is the characteristic Lan- 
dau level (LL) spacing depending on the magnetic length 
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Figure 8 Schematic of surface Landau level spectrum d52l ). 
For particle-hole symmetric case, only two dipole transi- 
tions N -> N + 1 (intraband) and —N -> N + 1 (in- 
terband) contribute to Hall conductivity, yielding the half- 
integer quantized a xy = (e 2 /h)(N + 1/2) d64]); N is the 
index of the highest occupied Landau level d59l ), Ep is the 
Fermi level. 



Ib = \/h/\eB\. The solutions for the LLs and eigen- 
spinors are given by 



e„ = sgn(n)h(2 B \/2\n\, \n) 



-iC^ n \\n\ - 1) 
C ln\\n\) 



(52) 

V2, (53) 




f-o = 0, |0) 



10} 



, C t0 = 0, Cm = 1, n = 



In the LL basis the Kubo conductivity tensor is 



ih y-v f n - f n > (n\j a \ri)(n'\jp\n 



2tt£ 2 b 



\ e„ - e n i huj + e n - e n > +ih/r' 

(54) 

where j a (/3) = ev& a (f}) is the a(/3)-component of the 
surface current operator (a, /3 = x, y), and /„ is the 
Fermi occupation number of the nth LL. Introducing oper- 
ators <r± = (a x ± ia y )/2 we find that the current matrix 
elements (n'\j x \n) = ev[(n'\a-\n) + (n'\a+\n)] and 
(n'\jy\n) = iev[{n'\a-\n) — (n'\a+\n)], obey the dipole 
selection rules 

(n'\a±\n) = ±iCi mn C^i) n '6\ n >^ n \ ±1 , (55) 

allowing only transitions between LLs with \n'\ — \n\ ± 1. 

4.2 Half-integer-quantized Hall conductivity cr xy . 
Role of particle-hole symmetry. We are interested in 
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the Hall component of the tensor, 



2 \ " 



(fn fn') 



\(n'\a + \n)\ 2 - \(n'\a.\n)\ 2 
(e n , - e n f/fi 2 + (r-i - iuSf 



in the QHE regime \/2Hbt = 4Ro<Ji 



(56) 



\B\ 



> 1. The 

realization of this regime depends on the classical surface 
conductivity a s = (e 2 /h)(kp£/2) and surface carrier den- 
sity n s (i?o = h/2e 2 and <Pq = h/\e\ are the resistance and 
magnetic flux quanta). At zero temperature LLs below the 
Fermi level Ep are occupied, 

-N h < n < N, f n = 1, (57) 

while those above Ep are empty (see also Fig.[8j: 




S(T) 



Figure 9 DC Hall conductivity a xy [given by Eq. d56l > for 
u = 0] in units of e 2 /h versus magnetic field B in Tesla. 
Inset: same dependence in high-field region; surface state 
parameters are Ro<7 s = 25, n s = 3 
330 meVnm, and temperature T = 2K. 



7V+l<n<iV p , f n =Q, 



(58) 



where N is the index of the highest occupied LL given by 



N = Int 



El 



\2h 2 0% 



Int 27m, 



Int 



n s $o 
\B\ 



(59) 

where Int(...) denotes the integer part. We also introduce 
the number of the hole LLs in the valence band, Nh, and 
the number of the particle LLs in the valence band, N p , 
both counted from the charge neutrality point. Using Eqs. 
d57l i and d58l we can write the Hall conductivity (T56b as 



N 



E E 



|<n'|<7 + |n>| 2 -|<nV_|n}| 5 



n'=N+l n=-N 

N p -N-l 



E E 

n'=N+l n=-N h 



( £n , _ en )2 /n 2 + (r _ X _ luj)2 

|Ha + H| 2 -|(^_H| 2 



(60) 



,(61) 



where indices n and n' run over occupied and empty LLs, 
respectively. In particular, the first sum in Eq. (l60b contains 
both /nfraband (within the conduction band) and /nferband 
(from the valence to the conduction band) transitions from 
the occupied LLs in the symmetric interval —N<n<N. 
The non-zero contributions come from the intra band tran- 
sition N —> N+ 1 with matrix element (N + 1 |cr+ \N) and 
from th e interba n d transi tion — N — > N + 1 with matrix el- 
ement (N + l\(J+ \ — N). The second sum in Eq. doTT i in- 
volves the interband transitions from the rest of the occu- 
pied (hole) LLs —Nh < n < —N — 1 to empty (particle) 
LLs N + 1 < n' < N p . This sum is nonzero only for 
the particle-hole asymmetric spectrum with different total 
numbers of the particle and hole LLs, N p ^ Nh- To see 



this we change n — > —n and then use |(n'|er_| — n}\ 
|(n'|(T_|n)| = | (n | a+ 1 n')|, which yields 



1/4 



(ejy+i - e N y/Ti z + (t- 1 - iujf 



1/4 



{e N+1 +e N y/h 2 + {T-^-iujy 

N p N h 



v v \(n'\a + \n)\ 2 -\(n\a + \n')\ 2 



(62) 



■(63) 



The summand in Eq. d63l is odd under exchange n <-> nf 
and, therefore, vanishes for N p = Nh. This means that the 
interband transitions induced by the spin-raising (er+) and 
-lowering (cr_) operators between the particle-hole sym- 
metric LLs n < — N — 1 and n' > N + 1 exactly com- 
pensate each other. In the QHE regime the dc conductivity 
a xy = a xy (ui = 0) d62l is 

^2 



1/4 



1/4 



(VN + T-VN) 2 (VTf + i + VN) 



N - 



N n = N h 



(64) 



Thus, the half-integer Hall conductivity d64l reflects the 
particle-hole symmetry in two ways. First, it involves the 
transitions from the symmetric particle and hole LLs, N 
and — N, to the lowest empty LL JV+1, and, second, inter- 
band transitions from deep LLs —Nh, —N— 1 to higher 
empty LLs N+l, N p do not affect <j xy provided that the 
conduction and valence bands contain equal LL numbers, 
N p = Nh- In Fig. [9] we plot the dc Hall conductivity <j xy 
in the whole magnetic field region and at finite temperature 
T, using Eq. dS&b with uj = 0. 

4.3 Particle-hole asymmetry N p ^ N h . Accu- 
racy of the half-integer quantization. It is instructive 
to discuss now the Hall conductivity a xy [see Eqs. d62l and 
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a)N p 



> N. 



b) N p<N h 




Figure 10 Particle-hole asymmetric Landau spectrum with 
(a) N p > N h and (b) N p < N h . Additional allowed inter- 
band transitions —Nh — > Nh + 1 (a) and 



-N p -1 



N n 



(b) give rise to non-universal corrections to the half-integer 
quantized Hall conductivity [see also Eqs. (|67| |, d68l l. and 



d63l ll for the particle-hole asymmetric LL spectrum with 
N p > Nh. Because of the odd summand in Eq. d63l ). the 
external sum involves the particle LLs Nh + N p : 



N ■ 



E 

%'=N h +ln= 



1 
2 

JV-1 



E 



+ 2(hf2 B ) 2 x 

\Wa + fa)\*-\Wa-W 



-N h 



y'- 



(65) 



(66) 



The last term contains an additional allowed transition 
from the deepest hole LL —Nh to the empty particle 
state Nh + 1 (see also Fig. [TOh ) with matrix element 



(N h + l\a + \ - N h ). Therefore, for N p > N h we have 



1/4 



(67) 



(VNh+T+VNh) 2 , 

In the similar way we obtain a xy for the particle-hole 
asymmetric LL spectrum with N p < Nh'- 



'.ry 



N 



1/4 



(68) 



In this case there is an additional allowed transition from 
the hole LL —N p — 1 to the highest empty particle state 
N p induced by the spin-lowering operator cr_ (see also Fig. 
[TOb). due to which the correction to the half-integer TV + 
1/2 is negative. To summarize, for large N p , Nh ^> 1 the 
Hall conductivity is given by 

, 2 



>xy 



e 

77 



N ■ 



sgn(JV p - N h ) 
l&-mm(Np,N h y 



(69) 




B(T) 



Figure 11 Two-surface DC Hall conductivity a xy [see 



Eqs. (l56l l and ( TTOb for ui = 0] in units of e 2 /h versus 
magnetic field B in Tesla. Inset: same dependence in high- 
field region; top (t) and bottom (b) surface state parameters 
are i? cr* = 25, n* = 3 • 10 11 cm- 2 and R cr b s = 12.5, 
n b s = 1.5-10 11 cm -2 . For both surfaces fiv = 330meV-nm 
and T = 2K. 



where correction 6 determines the accuracy of the half- 
integer quantization. 

4.4 Two-surface model. Since in a 3DTI two surface 
states, one at the top and one at the bottom of the sample, 
can contribute to transport II 101 . it is worthwhile to briefly 
discuss the behavior of the net Hall conductivity: 



(70) 



where the conductivities of the top and bottom surfaces, 
a x 'y, are both given by Eq. d56i l with, generally, differ- 
ent classical conductivities cr*' 6 and carrier densities n*' b . 
Such difference may arise when one surface faces a sub- 
strate, while the other - vaccuum II 101 . Although, sepa- 
rately, each surface exhibits the half-integer Hall quanti- 
zation, the plateaus in the net conductivity deviate from 
(e 2 /h)(N + 1/2) and are irregular, e.g. Fig. QT| shows 
an unusual plateau sequence e 2 /h,2e 2 /h,4e 2 /h,5e 2 /h... 
The QHE with unusual odd and even plateaus has been 
observed in the 3DTI HgTe where a bulk energy gap is 
induced by epitaxial strain 



1 1 1 1 . The bulk energy gap is 
a necessary prerequisite for the observation of the surface 
contribution in transport which is, otherwise, dominated by 
bulk carriers |97ll98ll99l[l00l . 

4.5 AC conductivities. Classical cyclotron reso- 
nance. Here we briefly discuss the AC Hall <r xy (uj) and 
diagonal a xx {ui) conductivities. At zero temperature T = 
for the particle-hole symmetric LL spectrum, they can be 
obtained from the Kubo formula ( T54b as 



1/4 



1/4 



(71) 



(^ / ^+ T +^) 2 + (7f^)' 
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e 2 1 — iujT 
h s/2£} b t 



(72) 



1/4 






( 1—iujT \ 


(y/N + 1 - y/N) 


(y/N + 1 - y/N) 2 + 
1/4 








f 1 — IUJT \ 


(y/N + 1 + y/N) 


(y/N + l + y/W) 2 + 





E 



1/2 



n>N+ 



1 s/n+l + 



In Eq. d72b the first and second terms arise from the intra- 
band (N -> N + 1) and interband (— JV ->• N + 1) transi- 
tions (see also Fig.O, while the third term accounts for the 
interband transitions from the hole LLs — iV— 2, — N— 1 
to the particle LLs N + 1, N + 2, .... The latter transitions 
contribute to a xx regardless of the presence (or absence) of 
the particle-hole symmetry. 

For weak magnetic fields when 



TV 



hn s 
~~eB 



> 1, 



(73) 



the main contribution to Eqs. (|7TT i and d72l comes from the 
intraband N — > N + 1 transition, leading to the classical 
AC conductivities 



(I2 c t) 2 + (1-zwt) 2 ' 

1 — ibJT 



(74) 



(75) 



(f? c r) 2 + (1 - jwt) 2 ' 

where fi c is the cyclotron frequency corresponding to a 
resonant transition between the highest occupied (N) and 
lowest unoccupied (N + 1) LLs (see also Fig. [8}: 



2eB 



(y/WTT-VN) 



eBv eBv 



(76) 



h hkp Ep 

Unlike the quadratic -dispersion case, fi c depends on the 



surface carrier density, f], 
wave- vector kp 



through the Fermi 



y/A-KUs, while v is the n s -independent 
band structure parameter. 

Complemented by the Maxwell equations, the magne- 
totransport theory discussed in this section lays the basis 
for magneto-optical spectroscopy of TIs, which is currently 
in the focus of both experimental lfTonil02U70lll 12111 181 
T221 and theoretical 1 1 2411 1 29lfT^[T^[n7lfB^fT62l[T74l 
research. In particular, the Faraday and Kerr effects have 
been observed in the TI materials fT02l7 WT2TITT8l[T22l . 
In high-mobility strained 3DTI HgTe, the Faraday effect 
reveals the QHE oscillations with low LL indices n ~ 1 
11221 . The broken time-reversal symmetry of the surface 



s 




s 


surface state 




surface state 










3DTI 
















E~ x 



Figure 12 (a) S/TI/S Josephson junction and (b) schematic 
energy diagram of the junction, showing Andreev bound 
states (ABSs) within induced superconducting gap Ai n on 
TI surface (see also text). 



QH states has been predicted to give rise to rich mag- 
netoelectric phenomena specific to axion electrodynamics 
0124M129l fT8l fT34lfT54l . 

5 Superconducting Klein tunneling on topologi- 
cal surfaces. Recently, experiments on the TIs have ad- 
vanced to the investigation of superconducting junctions 
involving 3DTIs as a weak link between conventional (sin- 
glet s-wave) superconductors (Ss) 11 13111 14|[1131fTT6lfTT7l 
11211 (see also Fig. [T2l . The novelty of such S/TI/S junc- 
tions is intimately related to the electron spin helicity: since 
for the helical states the spin-rotation symmetry is broken, 
a conventional singlet s-wave S (e.g. Nb or Al) is expected 
to induce not only the singlet s-wave pairing, but also the 
unconventional triplet p-wave pairing on the surface under- 
neath the superconductor Q123lll30IIT32lll49lll50lfT70l . 

The key question that arises at this stage is what are 
the potentially observable phenomena that could serve as 
a smoking gun for the p-wave superconductivity in S/TI/S 
systems? One of such phenomena, widely discussed in lit- 
erature, could be the formation of Majorana bound states, 
topological midgap states that are known to appear on 
edges (or in a vertex core) of a p-wave S (see e.g. Refs. 

The Majorana 

bound states are expected to give rise to unconventional 
Josephson effects 1 1831118411177112091 . dc superconduct- 
ing tunneling and current noise l40l . For the detection 
of the Majorana bound states it is essential, however, to 
avoid hybridization between different Majorana states. In 
Josephson junctions this requirement is not easy to fulfill 
because when two Ss are brought in contact their midgap 
edge states would normally hybridize to become a pair 
of Andreev bound states (ABSs) with finite energies de- 
pending on the phase difference between the Ss (see e.g. 
Ref. 1 184] for the case of intrinsic p- and d-wave Ss). In 
such a situation another closely related phenomenon - the 
superconducting Klein tunneling - can serve as a smoking 
gun for the p-wave superconductivity in S/TI/S junctions. 
Below we discuss this in more details. 



Copyright line will be provided by the publisher 



pss header will be provided by the publisher 



13 



5.1 Generalized Fu-Kane model of supercon- 
ducting proximity effect in 3DTIs. We begin by re- 
viewing the superconducting proximity effect for a single 
lateral contact between a conventional singlet s-wave S 
and the surface of a 3DTI (see also Fig.[T2l. Such a hybrid 
system was first considered by Fu and Kane II123II . In their 
approach the proximity effect on the surface is described 
by a singlet pairing potential treated phenomenologically 
as an energy-independent constant. On the other hand, mi- 
croscopic approaches (e.g. McMillan's model B2101 ) allow 
for a more general energy-dependent description of the 
proximity effect in terms of the Green's functions of the 
S. We will therefore follow McMillan's model 02 101 and 
its adaptations to low-dimensional systems (see e.g. 1121 II 
212 213 214|). In this model, the coupling between the 
systems is described by a tunneling Hamiltonian, allowing 
one to calculate the Green's function of the normal system, 
Gk, by summing up relevant Feynman diagrams generated 
by the tunneling Hamiltonian. The superconducting prox- 
imity is accounted for by a tunneling self-energy 27(e) in 
the equation for G^: 



[e-H k - 27(e)]G k = I, ff k = 



h k 

-h* 



(77) 



where is the Hamiltonian of the surface state in 2 x 2 
Nambu (particle-hole) representation, with = fiver ■ k — 
Ep being itself a 2 x 2 matrix in spin space. The self-energy 
27(e) is also a matrix in Nambu space with the following 
structure 112101: 



27(e) 



-£T(e)<7o A N (e)ia y c lx 
-A N (e)ia y e~ tx -ir(e)a 



(78) 



Its off-diagonal elements yield the induced singlet pairing 
potential (\ is the phase of the superconducting order pa- 
rameter), while the diagonal elements in Eq. ( 1781 account 
for the spectrum shift due to the tunneling: 

A N {e) = ir N f s (e) = ir N p4=F : (79) 



A 2 



r(e) = r N g s (e) = P N 



r=rrt 2 N s . (80) 



Here the tunneling energy scale is given by r N [it deter- 
mines the normal-state escape rate into the S, t is the tun- 
neling coupling strength, and N s is the normal-state den- 
sity of states at the Fermi level in S]. / s (e) and g s (e) 
are, respectively, the condensate and single-particle quasi- 
classical Green's functions of the S [A s is the gap energy 
in S]. From Eqs. (1771 and ( 1781 we explicitly find the effec- 
tive surface-state Hamiltonian with induced pairing: 



Hi 



eff 



27(e) 



(81) 



hvrr - k- E F -ir(e) A N (e)ia y e tx 
-A N (e)ia y e~ lx fiver* ■ k + Ep — i-T(e) 



Consequently, one can calculate all the Nambu matrix ele- 
ments of the surface-state Green's function: 



Gv 



Gn|k G 12 |k 

G21|k G*22|k 



(82) 



where Gii( 22 )ik and Gi 2 ( 2 i)|k are the particle (hole) and 
condensate Green's functions, respectively. We will briefly 
discuss the spin structure of the condensate function Gallic 
as it reveals induced mixed s + p-wave superconducting 
correlations: 



G 2 i 



1 

if) 



( a t-k a tk) 

K-k4k> 



2-(^0 



•kl 



(4-k a j k ) 

(4-k4k) 



iOyA, 



(6 + iT(e))s 



ft v 2 (k 



k F ) 2 -Al(e) 



(83) 



(84) 



Here we compare the general spin structure of G2i|k 
involving the ground-state expectation values of all time- 
ordered pairs of the creation operators a) with the explicit 
solution (l84l obtained from Eq.|77] Note that in addition to 
the singlet component ~ ia y , the spin-helicity cr ■ k gener- 
ates the spin-triplet p-wave component ~ cr ■ \do~ y with the 
same strength A N [k is the unit vector in the momentum 
direction on the Fermi surface]. Additionally, from the de- 
nominator of Eq. d84l we see that an isotropic energy gap 
Ai n is induced in the surface spectrum (see also Fig.[T2l: 



A; 



\A N (o)\ = r N 



e,T„ < A s 



(85) 



The s + p - wave proximity effect is robust against potential 
impurity scattering II 1491 . 

5.2 Superconducting Klein tunneling and topo- 
logical Andreev bound states. We consider now a 
short weak link with length L <C fiv/r N , choosing the 
superconductor phases on the left and right as \l =0 
and xr = X- At the simplest level the junction can be 
described by the following real-space equations: 

[H eff (r) + T z a U8(x)]V(r) = ei'(r), (86) 
&(x,y = 0)=V(x,y = W), (87) 

for the Nambu spinor SP = ^14.), (S^fi ^2i)] T which 

combines the particle (^1^, ^i^) T and hole (#2t) ^2i) T 
spinors. The transverse periodic boundary conditions ( ISTl i 
can be realized in a surface array of Josephson junc- 
tions. In order to model the interface scattering inside 
the weak link we introduce the potential U S(x) (see Ref. 
Ill 811 for the scattering-free case). The solution of Eq. 
(l86l can be sought as the sum of independent channels 
if(r) = Y.k n ^(xy^y/VW, with k n = 2-im/W and 
n G Z. Searching for ABSs with !f r fe n [x ±00) -4 Owe 
obtain the eingevalue equation: 




rri ■ 2 X 

T n sm -, 



(88) 
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-A„ 




Figure 13 Topological gapless ABSs E^(x) d90l 
They propagate perpendicular to junction barrier without 
backscattering. These states are orthogonal and 47r peri- 
odic (see also text). 



Figure 14 Nontopological ABSs E^(x) d93l correspond 
to oblique incidence at junction barrier. Scattering from 
barrier results in gapped 2ir -periodic spectrum (see also 
text). 



where T„ is the n's channel normal-state transparency: 



1 - {K/kpf 



i-(k n /k F ) 2 /(i + v?y 



U 
hv ' 



(89) 



with \k n \ < k F , determining the number of the open chan- 
nels in the junction, n c h- Eq. d88l accounts for the energy 
dependence of both A N (e) and -T(e). 

Superconducting Klein tunneling and topological ABSs 
occur in the n = channel, as it propagates perpen- 
dicularly to the junction barrier and is protected against 
backscattering (To = 1, see also Fig. [13). We obtain the 
ABS spectrum from Eq. (l88l using expansion in small pa- 
rameter 7 = r N /A s < 1: 



Et =0 (x) « ±r„ 



(l- 7 + 7 ')cos| + ^-cos a ^ 



.3 X 



(90) 



These states have orthogonal wave functions. Their spinor 
structure is simplest at zero energy: 



2^ ' 

^—ikpx—x I £ 

2V? ' 



(91) 
(92) 



where £ = hv/T N is the healing length. The orthogo- 
nality of d9"TT l and d92l prevents mixing of the counter- 
propagating modes, yielding the gapless 47r-periodic ABS 
spectrum ( f90b . Moreover, the particle d9"TT l and hole (l92l 
midgap states are constructed of the same two spinors 
[1, 1] T and [1, — 1] T describing two degenerate orthogonal 
Majorana modes. Similar midgap states were found in the 
core of a vertex on the 3DTI surface 01231 (see also reviews 
mi2l fT70lfT69l for more details). 

We should emphasize the difference between the topo- 
logical ABSs (|90T > and those appearing in p-wave Joseph- 
son junctions (see e.g. Ref. [184]). First, because of the 
superconducting Klein tunneling the ABSs d9"0l l are com- 
pletely independent of the details of the junction barrier. 
Second, the energy dependence of the proximity effect 



leads to the cubic and, in general, to higher odd-power 
terms in the spectrum. Additionally, due to the energy de- 
pendence of the proximity effect Eq. d88l has high-energy 
solutions appearing just below the superconductor gap A s 
(see thin blue curves in Figs. [13] and [141 these solutions 
will be discussed in detail elsewhere). 

The transport channels with k n ^ support non- 
topological ABSs because, at oblique incidence, scattering 
from the junction barrier generates an energy gap in the 
spectrum, making it 27r-periodic (see also Fig.fPfli: 



(1 - 7 + 7 2 )\/ 1 - T1 n sin2 (x/2)- 



+ ( 7 2 /2) [l-T„sin 2 ( X /2)] 



3/2 



(93) 



The topological d90l ) and usual d93l ABSs give rise to a 
peculiar AC Josephson current discussed below. 

5.3 Fractional AC Josephson effect. The AC 

Josephson effect is observed in voltage-biased junctions 
where the supercurrent oscillates with the frequency 
u> v = 2eV/h proportional to the bias voltage V and 
Cooper pair charge 2e 1121 5112 16l . In weak links between 
intrinsic p-wave superconductors, the An -periodic ABSs 
are expected to give rise to an unconventional AC Joseph- 
son effect at the fractional frequency lj v /2 (see e.g. Refs. 
11831118411 ). Leaving aside calculations (see Ref. ITMl for 
details) we present the result for the AC Josephson current 
in our proximity S/TI/S system: 

J(t) w J K-i) + Ji sin ^-t + J 2 sin ^-t + (94) 

Here the function Jo (lo v t) is the conventional AC Joseph- 
son current carried by the 27r-periodic ABSs ( |93l ), while 
the other two terms arise from the topological 47r-periodic 
ABSs ([90]l. Unlike the p-wave junctions 1183111841 , not 
only the main fractional frequency 0J v /2, but also Zuj v /2 
and, in general, higher odd harmonics appear due to the en- 
ergy dependence of the proximity-induced A N (e) d79l and 
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r(e) (f80b . The estimates of Jo, Ji and J 2 are given by 

Jo-^f^(l-7 + 7 2 ), (95) 

We can conclude that favourable conditions for the obser- 
vation of the topological ABSs exist in narrow weak links 
with a small number of open channels n c h ~ 1. 
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